A Method To Find Quantum Noiseless Subsystems 
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We develop a structure theory for decoherence-free subspaces and noiseless subsystems that applies 
to arbitrary (not necessarily unital) quantum operations. The theory can be alternatively phrased 
in terms of the superoperator perspective, or the algebraic noise commutant formalism. As an ap- 
plication, we propose a method for finding all such subspaces and subsystems for arbitrary quantum 
operations. We suggest that this work brings the fundamental passive technique for error correction 
in quantum computing an important step closer to practical realization. 
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Introduction. — The problem of controlling and main- 
taining properties of quantum systems which are in con- 
tact with an environment has received considerable re- 
cent attention. Primarily these investigations have been 
driven by the need to better understand the special fea- 
tures of evolving quantum systems that distinguish the 
quantum computing paradigm. Of central importance 
in the field of "quantum error correction" is the require- 
ment for techniques to avoid and overcome the degrad- 
ing effects of decoherence. Early work in quantum er- 
ror correction included the realization that many phys- 
ical error models contain symmetries induced by the 
system-environment interplay. This led to the discov- 
ery of "decoherence-free subspaces" (DFS) and "noise- 
less subsystems" (NS) Q, S II H S 11 ■ (On occasion 
we shall refer to both notions jointly as "NS" .) In these 
schemes, subspaces - and subsystems in the more ab- 
stract case - are identified within a system Hilbert space, 
with the property that all initial states encoded therein 
remain immune to the errors of a quantum operation of 
interest. Experimental efforts [a, 111 llfj, have affirmed 
the viability of this "passive quantum error correction" 
(PQEC) technique. 

It is also becoming clear that the NS formalism is appli- 
cable beyond the realm of quantum error correction. In 
quantum communication and cryptography, for instance, 
NS have been used as vehicles for avoiding noise (l2^ : 
this may lead to practical applications of NS in the near 
future. Further, NS are ideal for determining how to 
achieve distributed quantum information processing in 
the absence of shared reference frames . The NS con- 
cept has also arisen in recent analysis of black holes , 
and quantum gravity [J^, where NS are used to iden- 
tify the relational symmetry-invariant physical degrees 
of freedom in the quantum causal history framework. 

There is an obvious advantage to PQEC in the context 
of quantum computing. If a quantum operation (or chan- 
nel) is found to possess NS, then, by taking care at the 
initial encoding stage, the need for active error correc- 
tion after the fact is minimized. However, this protocol 



has a notable drawback. While substantial analysis has 
been carried out in important special cases, the protocol 
lacks a general method to find NS for arbitrary quantum 
operations. It is our belief that the PQEC approach will 
play a substantive role in quantum computing devices, 
and in applications beyond quantum computing, only if 
some sort of general approach for finding NS is derived. 

In this paper we propose such a method. Specifically, 
we develop a structure theory that shows precisely how 
properties of a quantum operation as a superoperator de- 
termine its NS structure. Moreover, if an operator-sum 
decomposition of Kraus (or "error") operators £ = {E a } 
for a channel £ is known, we show how algebraic proper- 
ties of the operators E a determine this structure. This in- 
formation naturally leads to the aforementioned method. 
Our analysis utilizes the framework for NS recently in- 
troduced under the umbrella of "operator quantum error 
correction" in [3113 (see also Q). 

As a consequence of this work, we suggest that the 
fundamental passive technique for error correction in 
quantum computing has been brought an important step 
closer to practical realization. Let us discuss these points 
further through a pair of illustrative examples, full details 
of the theory will be provided below. 

First we consider a simple example. Let H = C 2 <g> C 2 
be the combined system Hilbert space for two spin-^ par- 
ticles. Let {|00), |01), 1 10), |11)} be the associated basis. 
Consider the channel £ = {Z\, Z^\ where Z\ = Z ® la 
andZ 2 = t 2 ®Z, with the Pauli matrix Z = |0><0| — j 1><1 1 . 
Then the action of £ on a density matrix p on H is given 
by £{p) = \{Z x pZ\ + Z 2 pZ\). This channel has no non- 
trivial NS. The key point is that the "noise commutant" 
A' = {Zi,Z2}', which is the set of all operators on H 
that commute with both Z\ and Z2, only contains the di- 
agonal matrices with respect to the standard basis; i.e., 
the matrices corresponding to classical states. On the 
other hand, suppose our channel is T — {U±, U2}, where 
U k = UZ kl k = 1,2, and U = I4 — 2 jll)<ll|. In this 
case, the noise commutant A' — {U\, U2] 1 contains a sin- 
gle qubit NS (in fact it is a DFS). Indeed, any operator of 
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the form a = a|00)(00| + 6|00)(11| + c|00)(ll| + d|ll)(ll|, 
a, b,c,d £ C, belongs to A' and satisfies £{a) = cr. As dis- 
cussed below, that this NS is also fixed by £ follows from 
the fact that £ is unital, or bistochastic; i.e., £ (1) = 1. 

As a new example of NS, and one that will also be 
discussed further below, we consider an error model first 
discussed in |l9| in the context of active error correction. 
In this case the channel £ — {Eq, Ei,E%} acts on 2-qubit 
space and has three Kraus operators given by 

E a = a (|00)(00| + |ll)(ll|) + |01)(01| + |10)(10|, (1) 
E 1 = /3(|00)(00| + |10)(00| + |01)(11| + |11)(11|), (2) 
E 2 = /3(|00)(00|-|10)(00|-|01>(11| + |11)(11|), (3) 

where q is a scalar < q < 1 with a = y/1 — 2q and 
(3 = y/q/2. One can check that £{t) = Lo E i E \ ^ ^ 
and hence £ is non-unital. As we show below, the noise 
commutant here A 1 — {Eo, Ei, E%Y supports a single 
qubit NS that is not fixed by the action of £ . Further, 
there is another NS for the channel, in fact a DFS, that 
is not contained in the noise commutant. In particular, 
if we define the projector P = |01)(01| + |10)(10|, then all 
operators supported by P are fixed by £\ that is, £(a) — 
a for all a = PaP. However, these operators do not 
belong to the noise commutant. For instance, notice that 
EiP = 7^ PEi for i = 1,2. Hence, this error model has 
a NS inside its noise commutant that is not fixed, and a 
fixed DFS that is not contained in its noise commutant. 

Thus, one can ask, what is the underlying phenom- 
ena that produces noiseless subsystems? The previous 
example indicates that we must consider more than the 
noise commutant and fixed point set for the map. As it 
turns out, the structure theory we derive for NS can be 
phrased in terms of more general operator algebras ob- 
tained in the same spirit as the noise commutant, and, 
alternatively, in terms of modified fixed point sets for the 
map. Therefore, our approach has the advantage of ei- 
ther being set in an algebraic context, or strictly in terms 
of properties of the superoperator. 

The rest of the paper is organized as follows. We next 
recall the NS framework. We follow this by proving a 
theorem that yields the structure theory, and then show 
precisely how it may be used to find NS. Optimality of 
the method is then established, and this is followed with 
a conclusion on possible future work and limitations. 
Noiseless Subsystem Framework. - Given a quantum 
operation (or "channel" ) , represented by a completely 
positive, trace preserving superoperator £ : B(H) — > 
B(7i) on a (finite dimensional ) H ubert space Tt, the NS 

protocol QESEHUHlISIIl seeks subs y stems 
Ti B (with dim7i s > 1) of the full system Hilbert space 
H = (Tl A <g> TL B ) © K, such that 

Va A Va s , 3t a : £{a A ® a B ) = t a ® a B . (4) 

Here we have written o~ A (resp. a B ) for operators in 
B(H A ) (resp. B{H B )). In terms of partial traces, Eq. © 



can be equivalently phrased as, 

(Tr A o£)(a) = Tr A (cr), Vcr = a A % a B . (5) 

Thus, to be precise, B is said to encode a noiseless subsys- 
tem (or decoherence-free subspace in the case dimTt A = 
1) for £ : B(H) -► B(H) when Eq. ® is satisfied. 

The basic questions we address are the following: (1) 
Is there a structure theory for such subsystems? (2) If 
so, can it be applied to derive a canonical method to find 
such subsystems for arbitrary quantum operations? Our 
answer to the first question is yes, and for the second we 
make a proposal that lends itself to the possibility of a 
computational algorithm. 

Structure Theorem. — Let £ : B(H) — * B(H) be a quan- 
tum operation. We shall write £ — {E a } when an error 
model for £ is known; i.e., the operation elements E a 
determine £ through the familiar operator-sum represen- 
tation £(a) = Y, a E aO-El [HIHI. 

The (full) noise commutant A' for £ is the set of all 
operators in B{TL) that commute with the operators E a 
and E\. The f-algebra A generated by the E a is called 
the interaction algebra associated with £. In the unital 
case (£ (i) = 1) it is obvious that every a £ A' satisfies 
£ (a) = a, and, in fact, every operator that is fixed by 
£ belongs to A' jl^l • Of course, in the general case the 
operator £"(11) may not be so well behaved, and all that 
can be said for operators a £ A' is that they satisfy 
£ (cr) = a£(l) = £(i)a. This equation is suggestive of the 
more general phenomena that must be analyzed to obtain 
NS for arbitrary quantum operations. Given a projection 
P in B(H), we shall make the natural identification of the 
subalgebra PB(H)P of B(H) with the algebra B{PH). 

Theorem 1 Let £ — {E a } be a quantum operation on 
B{TL). Suppose P is a projection on TL such that 

£{P) = P£{P)P. (6) 

Then E a P = PE a P, Vo. Define 

A' P := {a £ B{PH) : [a, PE a P] = = [a, PEIP}}. 

and, 

Fix P (£) := {cr £ B(PH) : £{a) = a£(P) = £(P)a, 

£(a^a) = j£(P)a, £{aa^) = o-£(P)(T^}, 

Then Fixp(£ ) is a \-algebra inside B{PTL) that coincides 
with the algebra A' P ; that is, 

Fix P (£)=A , P . (7) 

Proof. Let P be a projection that satisfies Eq. ©. Then 

< P ± E a PElP ± < P ± £(P)P ± = Va. 

Hence P^E a P = 0, or equivalently E a P = PE a P, Va. 
Let E at p := PE a P = E a P, Va. It is clear that Fix P (£) 
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contains the commutant (taken inside B(PTC)) of the op- 
erators {E a P , E'l p }. Let a e A' P . We are required to 

show that a commutes with the operators E a ^ P and E p . 

The properties £{a) = a£(P) = £(P)a, a = PaP and 
Eq. (jSJ are seen through a calculation to imply that 

£(o*v) - a^£(P)a = Yy> E lj$W> E lj>\ ^ °- 

a 

(This inequality may be regarded as a generalization of 
the Schwarz inequality for completely positive maps from 
(Hill.) Thus, given £(a) = a£(P) = £(P)a, and so 
£ (crt) = £{aY = ^£{P) = £(P)ert, it follows that 

£ (crV) = a i a£(P) iff oE\ p = E\ p a, Va, 

£ (era 1 ) = o-cr^ (P) iff crP a , P = E a>P a, Va. 

This completes the proof. □ 
Observe that the maximally mixed state P = 11 triv- 
ially satisfies Eq. 10, and the algebra A% coincides with 
the full noise commutant {E a ,El}. However, as dis- 
cussed above, the operator 5(11) may not have many nice 
properties. In general there may be other projections P 
that support larger noiseless subsystems. 
Noiseless Subsystems. - Let P be a projection that 
satisfies Eq. (0. The structure theory for f-algebras 
yields a unitary U on PTi such that 

UA'prf = 0(l rofe ®M n J, (8) 

k 

for a unique (up to reordering) family of positive integers 
mk,rik > 1. We have used M nk to denote the operator 
algebra B(C nk ), represented as matrices with respect to 
some orthonormal basis. Note that the algebra A' P may 
be regarded as a subalgebra of B(TC) simply by taking 
a direct sum A' P O np of A' P together with the "zero 
algebra" O np of n P x n P matrices on P 7i, where n P — 
dimW - J2k m knk- 

The algebra structure Eq. (JSJ) induces a decomposition 
of the subspace PTi as 

PH = Q)(H Ak ®H Bk ), (9) 

k 

where dim(7i j4fc ) = rrik and dim(7i Bfc ) = n^. Observe 
that the positive operator £(P) belongs to the commu- 
tant inside B(PTL) of A' P by definition. As this com- 
mutant has structure Ap := A P = (B k (M mk ® ~R nk ), it 
follows that there are operators a^ S B(Tt Ak ) — M mk 
such that £ (P) = £) fc a k ® ll„ fc . 

Now let p — l Ak <g) p Bfc belong to the subalgebra 1 A * <g> 
S(W Bfc ) of ^ = ® k (l Ak ®B(H Bk )). Then we have 

£(l Ak ®p Bk ) = £{p) = p£{P) = <r fe ®p B *. (10) 

But Eq. (gj) holds if and only if it holds for a A = 1 A 
[HI . Therefore, it follows from Eq. {TUJI that each of 
the subsystems TL Bk is noiseless for £ and the following 
result is established. 



Theorem 2 Let £ be a quantum operation onB{7i). Let 
P be a projection on 7i that satisfies Eq. and let 
PT-L = © fc (^' 4fc ® ^ Bfc ) be th e decomposition of PH in- 
duced by the ^-algebra structure of A' P = Fixp(£ ). Then 
the subsystems Ti Bk , with dim7i Bfc > 1, are each noise- 
less subsystems for £ . 

In fact, it follows that if the input states are restricted 
to the subspace Tt Ak <£> Tt Bk , then the corresponding re- 
striction of £ satisfies £(Pk{ )Pk) = £k ® where Pk 
is the projection of Ti. onto TL Ak ® 7i Bfc , £k is a quantum 
operation on B(Tt Ak ), and ids fc is the identity channel 
on B(H Bk ). 

The NS structure for a number of unital channels 
have been analyzed in detail. An extensively studied 
class of channels arise from "collective noise", which 
has a number of physical interpretations (see |26l E?| 
for a detailed analysis of this and related NS struc- 
tures). We note a connection with [23 which includes 
a decomposition for collective noise channels of the form 
£ = J2k(Dk ® idfc)(PfcoPfc), where the Pt are projections 
associated with a decomposition of the system Hilbert 
space induced by underlying representation theory and 
the T>k are depolarizing channels (see Eq. (22) of 28]). 
Interestingly, this may now be seen as a special case of 
the general form derived here. 

Let us return to the non-unital example discussed in 
the Introduction. A computation shows in this case that 
the full noise commutant satisfies A\ — {Eq, Ei, E2}' = 
II2 8> M2, and thus supports a single qubit NS. Indeed, if 
o~ G A\ is written as H2 ® tro, Co € M2, with respect to 
this unitary equivalence, then a calculation shows that 

But recall that the projection P = |01)(01| + |10)(10| 
defines a DFS for £ ; specifically, £{a) — a, Ver = PaP. 
Now we can see precisely how this DFS arises. Namely, 
£{P) = P satisfies Eq. © and thus we find a single qubit 
NS for £, with \ip) := |01) and \(j>) := |10), given by 

A' P = spaniel, W(<Pl \<f>M, MM} = M 2 . 

Notice that A' P is not contained in A'%, and thus this 
DFS would not be detected through an analysis of the 
full noise commutant alone. Further, while A' t and A' P 
have the same "size" from an encoding viewpoint (i.e., a 
single qubit) (2^, it is perhaps more convenient to work 
with A' P = M2, as it can be more easily isolated within 
the full system Hilbert space. 

In fact, this example gives an indication as to how ac- 
tive and passive techniques for quantum error correction 
can be combined to combat noise. Indeed, we have just 
noted that the subspace ||01), |10)| determines a DFS 
for £. On the other hand, in |l9j it was shown that ac- 
tive error correction may be used to overcome corruption 
by £ of the code subspace {|00), |11)}. 
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Optimality of the Method. — The previous two sections 
yield a canonical method to compute noiseless subsys- 
tems for a given quantum operation £ which can be suc- 
cinctly stated as follows: 



channels, etc. It is important to note that the full noise 
commutant captures all NS in this case. In particular, 
this means algebras A' P may not be contained inside A\ 
only in the non-unital case. 



(i) Compute the projections P such that Eq. © holds. Corollary 1 Let £ be a unital quantum operation on 



(ii) Compute the structure of the algebras A' P = 
Fixp(£) as in Eq. ©• 

Then, in the notation above, the subspaces TC Bk , with 
dimH Bk > 1, encode noiseless subsystems for £ via the 
operator algebras l Ak <g> B{H Bk ) 

A crucial final step in the process is to determine if 
this scheme captures all noiseless subsystems for £ . We 
next show that this is indeed the case. 

Theorem 3 Let £ be a quantum operation on B{TL). 
Suppose that H = (H A ®U B ) © K and that H B is a 
noiseless subsystem for £ as in Eq. p|). Let P be the 
projection of H onto H A <g> H B . Then £{P) = P£(P)P 
and the algebra A' P contains 11 <g) B(H. B ) as a simple, 
unital ^-subalgebra. 

Proof. Let be an orthonormal basis for TL A , and 

let {Pki — \ctk){ai\ ® 11 s } be the corresponding matrix 
units inside B{H A ) <g> 1 B . It was proved in [3, 113 that 
Ti B is noiseless for £ = {E a } as in Eq. 0} precisely when 
E a P = PE a P and there are scalars {\ a ki} such that 



PhhEaPl 



ll 



KkiPki Va,fc,Z. 



(11) 



Note that the projection P is given by P = J2k^ k ' 
where we have written P^ for Pkk- Thus we have 

£{P) = J2 E aP k El= ]T PtE^ElPv 



a , k 



a,k,l,l' 



— ^alk^al'kPlkPkl' — 



a,k,l,l' 



i.k,l,l' 



^alk^al'kPll' ■ 



Let a = IT* 1 <g> a B e IT <g> B(H B ). Then since the P kl 
commute with a — PaP we have 

£ (a) = £{PaP) = P£(PaP)P 

= J2 PkEaPvaPvElPi 

a.k,k' ,1,1' 

= ^akk' ^all' Pkk' & Pi' I 

a.k,k' ,1,1' 

= a£{P) = £{P)a. 

In particular, this implies (with a B = I s ) that £{P) = 
P£{P)P and that the algebra 1 A <g> B{H B ) is contained 
in A' P . It is clear that A' P and 1 H A ^)B(H B ) have the same 
unit P, and that 11 B(Ti B ) is a simple (i.e., contains 
no non-trivial ideals) f-subalgebra of A' P . □ 

We finish with a consequence for the unital case. The 
class of unital channels includes numerous physical er- 
ror models such as collective noise, randomized unitary 



B(H). If % = 11 ®B{H B ) is the algebra determined 
by a noiseless subsystem for £ as in Eq. then 21 is a 
subalgebra of the full noise commutant A\. 

Proof. Let P be the projection of Ti. onto TC A ® Tt B . By 
Eq. Q, there is a r A such that £{P) = £{i A <g> 11 s ) = 
t a ® 11 s . Since £ is a unital completely positive map, 
we know that t a is a contraction operator, and hence 
£JP) < P. Then in fact £{P) = P by Lemma 2.3 from 
[22J. Thus, it follows from Theorem|21 and the definition 
of A' P , that 21 C A' P is a subalgebra of Fix (£). □ 
Conclusion. — We have derived a structure theory for 
decoherence-free subspaces and noiseless subsystems that 
applies to arbitrary quantum operations. As an applica- 
tion, we have proposed a method to compute NS for any 
given operation. We expect that the method could be for- 
malized into a computational algorithm, as suggested by 
recent literature |3(i \'M\ \ which includes algorithms writ- 
ten to calculate operator algebra structures, but there 
are still details to work through. We plan to undertake 
this investigation elsewhere. 

We discussed a non-unital example in which the max- 
imally mixed state and a smaller projection support dif- 
ferent single qubit noiseless subsystems. We suggest that 
this work motivates reconsideration of the quantum chan- 
nels that appear in the literature, for the possible exis- 
tence of noiseless subsystems. We wonder about possible 
experimental implications of this work. It would also be 
interesting to investigate connections with other recent 
noiseless subsystem related efforts such as EE • 
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